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We have studied the effect of polyethylene glycol (PEG) on the 
liquid-liquid phase separation (LLPS) of aqueous solutions of 
bovine yD-crystallin (yD), a protein in the eye lens. We observe that 
the phase separation temperature increases with both PEG con- 
centration and PEG molecular weight. PEG partitioning, which is 
the difference between the PEG concentration in the two coexist- 
ing phases, has been measured experimentally and observed to 
increase with PEG molecular weight. The measurements of both 
LLPS temperature and PEG partitioning in the ternary yD-PEG- 
water systems are used to successfully predict the location of the 
liquid-liquid phase boundary of the binary yD- water system. We 
show that our LLPS measurements can be also used to estimate the 
protein solubility as a function of the concentration of crystallizing 
agents. Moreover, the slope of the tie-lines and the dependence of 
LLPS temperature on polymer concentration provide a powerful 
and sensitive check of the validity of excluded volume models. 
Finally, we show that the increase of the LLPS temperature with 
PEG concentration is due to attractive protein-protein interactions. 
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Polyethylene glycol (PEG) is a hydrophilic nonipnic polymer 
used in many biochemical and industrial applications, frue to 
its nontoxic character, this chemical can be found in cosmetics, 
food, and pharmaceutical products. The mild action of PEG on 
the biological activity of cell components explains the success of 
this polymer in biotechnological applications. PEG is commonly 
used for liquid-liquid partitioning and precipitation of biomac- 
romolecules (1,2). In protein crystallography, PEG is considered 
the most successful precipitating agent for the production of 
protein crystals, the crucial step for the determination of the 
molecular structure of a protein. All these applications make 
PEG by far the most widely used polymer in aqueous solutions 
of biological molecules (1, 3). 

Due to the extensive practical use of PEG as a precipitating 
agent for proteins, it is of fundamental importance to understand 
protein-protein and protein-PEG interactions in protein-PEG- 
water ternary systems. These interactions are often described in 
terms of a depletion force, which arises because the polymer is 
depleted in the region between adjacent proteins (4, 5). Deple- 
tion force models have been successful in describing the effect of 
nonadsorbing polymers on colloidal suspensions (6-10). 

Protein-PEG-water solutions have been investigated by sev- 
eral techniques (4, 5, 11-14). Small-angle x-ray- (4) and light- 
scattering (5) measurements generally confirmed that total 
protein-protein interactions can be described in terms of deple- 
tion effects. The microscopic interpretation used in the above 
studies is derived from colloid-polymer models. To model the 
effect of PEG, a depletion component is added to the original 
protein-protein pair interaction potential. This microscopic 
modeling, however, does not provide information about the 
actual protein-PEG interactions in the ternary system. Dialysis 
experiments, which describe the effect of PEG partitioning (i.e., 
the difference of PEG concentration) between two solutions 
separated by a membrane not permeable to protein molecules, 
give information on protein-PEG interactions only when pro- 



tein-protein interactions are absent (12, 15). Studies of the 
solubility of protein crystals in the presence of PEG have been 
carried out to investigate protein-protein and protein-PEG 
interactions. However, to interpret the results, assumptions 
about the properties of the crystal are required. Moreover, 
simple excluded volume models fail to describe protein solubility 
as a function of PEG molecular weight (16). The discrepancy 
between model and experiment becomes particularly large for 
PEG of lower molecular weights (11). Furthermore, the mag- 
nitude of the protein-PEG interactions changes from protein to 
protein in a way that does not correlate with the simple excluded 
volume scenario (11, 15). 

Liquid-liquid phase separation (LLPS) can be used as a direct 
tool for probing both protein-protein and protein-PEG inter- 
actions. It has been shown that the position of the liquid-liquid 
boundaries provides insight about the nature and magnitude of 
the protein-protein interactions (IT^Ci)." Ttfe ! ;a'dva(ti^g9 : of 
LLPS analysis with respect to measurements of protein solubility 
is that data interpretation requires no assumptions regarding 
crystal phase properties (11, 21); LLPS is important' also for trie 
analysis of nonequilibrium properties of the system. Because 
the nuc.ieation of protein crystals is erjtharicdd in proximity bf¥he : 
LLPS, this phase boundary is believed to be'ih^^licatecl hi protein 
crystallization ? kinetics (22; 23). ■• niiMWi.-'" ^r-:.:i>;! vl - 

We present an experimental investigation of the- effect' of low 
molecular weight PEG on the LLPS of aqueous solutions of 
bovine yD-crystallin (yD). yD is a member 1 of the 1 y^crystallin 
family of lens proteins, which are involved in cataract formation 
(24). The phase diagram of aqueous y-crystallins has been well 
characterized, and several members of the family, including yD, 
have attractive interactions leading to LLPS of the correspond- 
ing aqueous solutions (17, 20, 25, 26). yD is a protein with a 
simple globular structure (27) and is weakly charged at the pH 
value (^7) used in our measurements (28). The experiments 
were performed at a PEG concentration of «*5% in weight, 
which falls within the range of concentrations used for protein 
crystal growth (3). 

The LLPS properties of the protein-PEG-water ternary sys- 
tem are described by a coexistence surface in the phase diagram. 
This coexistence surface represents the LLPS temperature as a 
function of protein concentration, c\, and PEG concentration, 
C2. At a fixed temperature, the LLPS properties of the system are 
described by an isothermal coexistence curve, which gives the 
concentrations (cj, c\) and (c", c") of the coexisting phases, I 
and II. The partitioning of the components in the two coexisting 
phases is described by tie-lines connecting the points (c{, c\) and 
(c", c l 2 ) of the coexistence curve. The critical point, (cj, c£), is 
defined as the point of the coexistence curve where the condition 
(cj, c\) = (c", c") occurs. The location of this point as a function 
of the temperature is described by a critical line on the coexist- 
ence surface. 



Abbreviations: PEG, polyethytene glycol; LLPS, liquid-liquid phase separation; yD, 
yD-crystallin. 
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Materials and Methods 

Bovine -yD, which has a molecular mass of 20,700 g/mol, was 
isolated from 1- to 6-week-old calf lenses, obtained by overnight 
express from Antech (Tyler, TX). Pure yD was isolated and 
purified from these lenses by standard procedures described 
elsewhere (25, 29). The purity of the native sample was at least 
95%, based on both ion-exchange and size-exclusion HPLC. The 
purified yD was dialyzed exhaustively into sodium phosphate 
buffer (0.1 M, pH 7.1) that contained sodium azide (0.02%). 

PEG with molecular masses ranging from 200 to 3,350 g/mol 
(PEG200-PEG3350) and tetraethylene glycol (TEG, molecular 
mass 194.2 g/mol) were purchased from Sigma- Aldrich and 
used without further purification. Sample molecular mass and 
polydispersity were analyzed by electrospray ionization mass 
spectrometry at the Biopolymers Laboratory at the Center for 
Cancer Research at the Massachusetts Institute of Technology. 
Estimates of the average molecular mass were in agreement with 
the nominal values provided by the company. All polymers had 
a narrow molecular mass distribution with a polydispersity of 
about 1.05. 

Protein-PEG aqueous solutions were prepared as follows. 
Solutions containing dilute yD in phosphate buffer were con- 
centrated by ultrafiltration. When the desired protein concen- 
tration was reached, a known weight of PEG was added to the 
protein solution. The concentration of yD in the samples was 
determined by UV absorption at 280 nm, using the extinction 
coefficient value of 2.11 mg'^mbcm^ 1 (26). The concentration 
of PEG in the samples was calculated by using the mass of PEG 
and the total volume of the solution. The volumetric contribution 
of each component is obtained by multiplying the mass of the 
component by the corresponding specific volume, i.e., 0.71 ml/g 
for yD (30), 0.84 ml/g for PEG (15) and 0.992 ml/g for the 
buffer. 

The LLPS temperature for a given protein-PEG aqueous 
solution was determined by gradually lowering the temperature 
of the sample until clouding was observed. The cloud point was 
determined by the examination of the transmitted intensity- 
temperature profile as previously described by Liu et al. (31). The 
coexisting phases were obtained by quenching the sample at 
fixed temperature beneath the coexistence surface. If after about 
24 hours the two coexisting phases had not separated by gravity, 
centrifugation was used for the separation. The protein concen- 
tration in each phase was determined by UV absorption. To 
determine the PEG concentration in each of the coexisting 
phases, an aliquot of known weight was taken from each phase, 
and the PEG was separated from the protein by ultrafiltration 
(Amicon Microcon YM-10). The protein-free solution was 
weighed and the corresponding concentration of PEG deter- 
mined by using a standardized refractive index detector (Perkin- 
Elmer LC-30 RI). The concentration of PEG in each aliquot was 
then calculated from its measured value in the protein-free 
solution. To ensure that possible differences in phosphate con- 
centration in the two phases did not affect the results, the 
refractive index measurements were performed after isocratic 
elution of the protein-free solution on a size-exclusion HPLC 
column, by using sodium phosphate buffer (0.1 M, pH 7.1) with 
a flow rate of 1 ml/min. An Ohpak SB-802.5 HQ column (size 
8 mm X 300 ml) from Phenomenex was used for the PEG200 and 
PEG400 solutions. For PEG1000, the Superdex 75 HR 10/30 
column (size 10 X 300 mm) from Amersham Pharmacia Biotech 
was preferred for its higher resolution. The procedure was 
verified with protein-PEG aqueous solutions of known compo- 
sition. In all cases, the measured protein and PEG concentra- 
tions in the two coexisting phases were consistent with the 
protein and PEG concentrations in the original homogeneous 
samples. 
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Fig. 1. LLPS temperature at constant PEG concentration of ~50 mg/ml for 
the yD-PEG-water ternary system. The average molecular masses are: 200 
g/mol (■), 400 g/mol (•), 1.000 g/mol (♦), 1,450 g/mol (A), and 3.350g/mol 
(t). The solid curves are guides for the eye. The values for the yD-tetra ethylene 
glycol-water ternary system are represented by open squares (□). We draw the 
coexistence curve for the yD-water binary system (dashed curve). The vertical 
bars (j) locate the critical point. 



Results 

We describe the composition of the yD-PEG-water systems by 
the protein (component 1) concentration c\ and the PEG 
(component 2) concentration inside the protein-free volume, c^. 
If v is the protein specific volume and <f> = C\v is the protein 
volume fraction, then c-^ is equal to c 2 /(l - <£). The choice of 
C2s instead of c 2 to present our results is justified by the following 
argument. We consider as a reference system one in which the 
particles of component 2 have a negligible size and do not 
introduce specific chemical interactions in the system. In this 
reference case, the different values of c 2 between the two 
coexisting phases are entirely determined by the exclusion of 
component 2 from the volume occupied by component 1. Here, 
the condition of chemical equilibrium is expressed by the rela- 
tion: = Thus, in our systems, any observed difference in 
C2s between two coexisting phases is related to the finite size and 
specific chemical properties of the PEG molecules. 

In Fig. 1, we present our measurements of the LLPS temper- 
ature, T ph , for yD-PEG-water ternary solutions at constant PEG 
concentration of C2s 50 mg/mi and with several polymer 
average molecular masses: 200 g/mol (filled squares), 400 g/mol 
(circles), 1,000 g/mol (diamonds), 1,450 g/mol (upright trian- 
gles), and 3,350 g/mol (inverted triangles). We also report 
measurements of temperature of phase separation for the yD- 
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Table 1. The values of the parameters characterizing the ternary coexistence surfaces 





C2ss mg/ml 


To K 


OTph/ac,)^ 
10- 3 K mg-'-ml 


10- 3 -mg- 1 -ml 


(dTph/dCzs)^. 
K-^mg-'-ml 


Tl K 


PEG200 


54 


280.6 


3 ± 2 


-0.5 ± 0.1 


0.11 ±0.07 


275 ±4 


PEG400 


52 


286.5 


9 ± 1 


-1.1 ±0.1 


0.17 ±0.03 


278 ±2 


PEG 1000 


50 


299.3 


33 ±3 


-2.0 ± 0.3 


0.33 ± 6.05 


283 ± 3 



tetraethylene glycol-water at the same value of C2s (open 
squares). TEG essentially has the same molecular weight as 
PEG200, but, in contrast to PEG, is monodisperse. We see from 
the figure that the two sets of measurements, for TEG and for 
PEG200, overlap within the experimental error and conclude 
that the effects due to polydispersity are small. In Fig. 1, we also 
show the coexistence curve for the binary aqueous protein 
system (dashed curve) (29). The phase boundaries of the ternary 
systems are located at higher T ph than that of the binary system, 
and the difference increases with PEG molecular weight. The 
critical points, shown by vertical bars in Fig. 1, were identified by 
the examination of the transmitted intensity-temperature profile 
as previously described by Liu et al (31). In all cases, the value 
of protein critical concentration (about 290 mg/ml) is not 
affected by PEG within the experimental error. On the other 
hand, the slope of T P h at the critical point, (dTp^/dCi)*^, is 
positive and increases with the PEG molecular weight, in 
contrast with the zero slope required for the protein-water 
binary system. In Table 1, we report the values of the PEG con- 
centration, C2s, the critical temperature, T C7 and the slope 

(ar ph /dc0^. 

In Fig. 2, we present our experimental data (filled circles) of 
the coexisting values (c\ , c^) and (c?, c£) for the yD-PEG-water 
solutions having PEG average molecular weight: 200 g/mol 
(279.2 K, case A), 400 g/mol (281.7 K, case B), and 1000 g/mol 
(290.3 K, case C). The pairs of points representing the coexisting 
phases are connected by straight lines (tie- lines). The dashed 
curves are eye guides representing the coexisting curves. In all 
three cases, the coexisting phases were obtained by quenching 
samples with the same initial protein concentration of 300 
mg/ml but with different initial PEG concentrations. The 
temperatures for the measurements of the coexistence curves 
were chosen so that they lie above the critical temperature of the 
protein-water binary system but below the critical temperature 
of the corresponding ternary system in Fig. 1. Because the 
protein critical concentration, c\, is equal to the average con- 
centration (c\ + c l i)/2 in the limit of c\ — c? = 0, we determine 
c\ by plotting (c\ + c?)/2 as a function of cf - c?. The resulting 
values of the critical concentrations (vertical bars in Fig. 2) are, 
within the experimental error, the same as the value reported for 
the binary protein-water system. An interesting feature of the 
phase diagram is that the difference of PEG concentration 
between the two coexisting phases increases as the PEG molec- 
ular weight increases. To describe the effect of PEG molecular 
weight on PEG partitioning, we consider the normalized slope 
of the tie-lines at the critical point: (dc^/dc i)r ph / c ly This 
quantity can be obtained by plotting the incremental ratio (In 
— In C2s)/(c\ — c l i) as a function of c\ — c" and considering the 
limit of c\ — c° = 0. In Table 1, we can see that (dC2j/dci)^ ph /c|, 
increases with the PEG molecular weight. The protein volume 
fraction in the protein-rich phase is comparable to the protein 
volume fraction in the crystal, which is «*0.4-0.6 (ci «* 600-800 
mg/ml) (27, 32). The nonnegligible values of PEG concentration 
of the protein-rich phases reported for the yD-PEGlOOO-water 
system suggest that the crystal phase could accommodate poly- 
mer coils with molecular weight equal to or lower than 1,000 
g/mol. 



The behavior of the coexistence surface is described by the 
three slopes (dc^,/ dci) Tpb , (dT ph / dC2*) Cl , and (dTpt/dc i) C2M y 
which are related to each other by the mathematical relationship: 

It is important to observe that Eq. 1 is true only for ternary 
systems. Due to the presence of the two buffer salts (dibasic and 
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Fig. 2. Coexisting surfaces at constant temperature for the yD-PEG-water 
ternary system with PEG average molecular mass: 200 g/mol (279.2 K, case A), 
400 g/mol (281.7 K, case B), and 1,000 g/mol (290.3 K, case Q. The pairs of 
points representing the coexisting phases (•) are connected by the tie-lines 
(solid lines). The dashed curves are guides for the eye, and the vertical bars (]) 
locate the critical point. 
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monobasic sodium phosphate), our protein aqueous solutions 
contain more than three components. We know that salt parti- 
tioning is unimportant if PEG is not in the solution, because the 
protein-buffer system is well described as a binary system [for 
example, the critical point is at the maximum of the LLPS curve 
(26)]. Therefore, in the protein-PEG- buffer system, it is only 
because of the presence of PEG that salt partitioning may occur. 
We will show that Eq. 1 successfully describes our experimental 
results, which suggests that salt partitioning is insignificant 
insofar as the validity of Eq. 1 is concerned. 

We apply Eq. 1 on the critical line of the coexistence surfaces. 
Because the critical temperature, T c , of the -yD-PEG-water 
systems (with czs ~ 50 mg/ml) is higher than the critical 
temperature, 7^!, of the yD-water system (Fig. 1), the quantity 
(dTph/ dczOc! xs positive. At the critical point, the slope 
(dC2s/dCi)r ph of the coexistence surface is negative because it 
coincides with the slope of the tie-lines (Fig. 2). Eq. 1 shows that 
the slope (&T P h/ dc{) C2t is positive, as found experimentally (Fig. 
1 and Table 1). 

The slope (dT ph / dc^)^ can be calculated, by using Eq. 1, from 
the values of czy, (dTph/dc^l^, and (dc^/dc O^AL reported 
in Table 1. We have assumed that the quantity, 
(dc^/ dc i)r pb / c 2s, which is obtained from our measurements of 
coexistence curves, is not a function of the temperature and PEG 
concentration within our experimental domain. The calculated 
values of (dT ph / dC2s)l v reported in Table 1, are consistent with 
the experimentally observed dependence of the critical temper- 
ature as a function of PEG molecular weight (Fig. 1). Indeed, we 
can quantitatively compare the experimental change of critical 
temperature between the ternary and the binary systems, shown 
in Fig. 1 for each PEG molecular weight, with the change 
predicted from the calculated values of (dT ph / dC2s)t r Spe- 
cifically, we may estimate the critical temperature, Xj!, of 
the yD-water system from the expression: 7^ — T c — 
(dT P h/ dc^cxCzs- Thereby, we obtain the calculated values of 7^, 
which are reported in last column of Table 1. These are in 
acceptable agreement with the experimental value of 278.4 K 
(29). Hence, using only the information obtained from the 
properties of the ternary coexistence surfaces, we are able to 
estimate reliably the location of the binary liquid-liquid phase 
boundary. This information can be most valuable when LLPS of 
protein-water systems is not experimentally observable because 
it is inaccessible due to either the freezing of the system or the 
denaturation of the protein at high temperature, which is 
important for proteins with lower consolute critical points (33). 

Discussion 

Our previous analysis shows that the behavior of the coexistence 
surfaces at the critical point can be described by two distinct 
features: the slope of the tie- lines and the dependence of LLPS 
temperature on PEG concentration. To analyze the physical 
factors that determine these two features of the coexistence 
surfaces, we introduce the free energy of the system. We define 
the quantity F, representing the difference, at constant volume, 
K, and temperature, T, between the Helmholtz free energy of the 
virtually incompressible protein-PEG-water system and the 
pure water system. The changes of F due to the replacement (at 
constant volume) of 71 water moles by 1 mole of protein and y 2 
water moles by 1 mole of PEG are respectively described by the 
differences of chemical potentials, /xi = /utp rotc in — YiAWater and 
/X2 s mpeg — 72£Wer- A ternary incompressible system may be 
equivalently treated as a binary compressible system where 
and fia are the chemical potentials of the two effective compo- 
nents (17). The quantities yu\ and fi 2 will be indicated as the 
protein and PEG effective chemical potentials. 

It is convenient to introduce the reduced free energy/ (F — 
F°)/RTV, where F° is the standard free energy, and R is the 
ideal gas constant. If the concentration of PEG, C2, is relatively 



small, the reduced free energy is, to first-order approximation 
with respect to cz, given by the following equation: 

fc u c 2 , T) =/( Cl , 0, T) + c 2 ln(c 2 A0 + c 2 £c lt T). [2] 

In Eq. 2, f(c\, 0, T) is the reduced free energy of the binary 
protein-water system, whereas c% \n(c 2 /e) andc2^(ci, T) are the 
contributions to the reduced free energy associated with the 
replacement of water molecules by PEG molecules at constant 
V and T. The contribution of PEG to the ideal mixing entropy 
of the system is represented by the term C2 ln(c2/e), whereas the 
quantity c 2 £(ci, T) is the first term in a series expansion 
describing the change of the reduced free energy over and above 
the ideal mixing entropy. This approximation is acceptable 
because our experimental PEG concentrations are generally well 
below the threshold for polymer-polymer interpenetration, i.e., 
we are in the dilute regime (5). If we differentiate Eq. 2 with 
respect to the concentrations, c\ and C2, at constant Vznd T, we 
obtain the following expressions for the reduced effective chem- 
ical potentials: 

Ai(ci, c 2 , T) = Ai(ci, T) + Cz&e/dcJn . [3a] 

£ 2 (ci, c 2 , T) = In c 2 + £c u T), [3b] 

where & = (/x, - n?)/RT (with i = 1, 2), im(T) = 
(dF°/dCi) T ,v/Vy and the quantity £i(ci, T) = £i(ci, 0, T) is the 
protein effective chemical potential corresponding to 
the protein-water binary system. The two crossderivatives of the 
chemical potentials* (d£i/«te2)c„T and (dfa/ dc\) c ^ T i are both 
equal to (d£/dci)r- This quantity characterizes the protein-PEG 
interactions. Eq. 3a and 3b fully characterize all the experimental 
observable features of the coexistence surfaces. Specifically, they 
enable us to analyze the quantities (dcz/dci)^ the slope 
of the tie-lines) and (dT ph /dc 2 ) Cl . 

The Slope of the Tie-Lines. If we differentiate Eq. 3a with respect 
to c 2 at constant Ci and T, and Eq. 3b with respect to c 1 at 
constant fa and T, we obtain: 

(*h) =(11) = _ 2J*e£\ r41 

Experimentally, the quantity (Bc2/dC\)^ T can be found from 
the limiting value of the slope of the tie-lines at the critical point. 
If we apply Eq. 4 to our experimental tie-line, we observe that 
(djli/dc 2 ) CliT = (d£/dci) T is positive and increases with the 
PEG molecular weight. 

The slope of the tie-lines can be used to predict the effect of 
PEG on protein solubility. If protein crystals (solid phase) are in 
thermodynamic equilibrium with the liquid phase, pL X must be 
equal to its value in the solid phase, /if. For relatively low protein 
concentrations, fi\(c\ t T) In c\. Thus, using the equality 
of chemical potentials and Eq. 3a, we find: In C\ •** /if — 
(d€/dCi) T c 2 (11). This relationship gives the protein solubility, 
cu as a function of PEG concentration, c 2 . The dependence of 
/li on Ci and c 2 is expected to be very weak. Thus the fractional 
decrease of protein solubility with PEG concentration is fully 
determined by the value of (d£/dci)r, which is proportional to 
the slope of the tie-line (see Eq. 4). Therefore the magnitude of 
(dc 2 /dCi)^ jT /c 2 characterizes the effectiveness of a protein 
precipitating agent such as PEG. 

The Excluded Volume Model of Protein-PEG Interactions. We now 

analyze the experimental coexisting values (</>*, c\) and (<f> n , c\ l ) 
to gain insight about the physical factors that determine the 
quantity £ and its dependence on C\. We apply a simple excluded 
volume model that directly relates PEG partitioning to the 
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fig. 3. Dependence of the size ratio, q, as a function of the square root of the 
PEG average molecular weight, M n . The filled circles Ware the value obtained 
by fitting the experimental coexisting compositions with Eq. 2; the open circles 
(O) are the values calculated from the polymer gyration radii and the protein 
molecular volume. The dashed lines are guides for the eye. 

difference in free volume fractions between the two coexisting 
phases. If we assume that PEG molecules can be described as 
ideal polymer coils and that protein molecules have a spherical 
shape, then due to steric hindrance, each protein will be sur- 
rounded by an adjacent region where the centers of mass of the 
coils are excluded, and the width, S, of the resulting depletion 
layer will be proportional to the gyration radius, R g , of the 
polymer coil (4-14). If a is the volume fraction available to the 
centers of mass of the coils, the condition of chemical equilib- 
rium £2(^1, c\, T) = £2(^1, C2, T) becomes c\/ol 1 = c^/a 11 , 
where C2/0C is the polymer concentration in the free volume 
(9, 34) and the free volume fraction, a, is equal to exp(— f) (see 
Eq. 3b). 

It has been theoretically shown that, in the case of ideal 
polymer coils and relatively large hard spheres, 5 = l.lR g (35). 
An approximate expression for the free volume fraction, a, as a 
function of the hard sphere volume fraction, </>, is given by the 
well established scaled particle theory (34, 36, 37): 

a = (1 - 4>)exp(-^7, - B V 2 - Ct, 3 ), [5] 

where tj = <f>/(l - <f>),A = 3q + 3q 2 + q 3 ,B = 9q 2 /2 + 3q 3 , 
C = 3q 3 9 and q = 8/R is the depletion layer thickness 
normalized with respect to the hard sphere radius, R. The 
prefactor, 1 — </>, in Eq. 5, is the volume fraction not occupied 
by the spheres, whereas the exponential factor describes the 
effect of the depletion layers. When <#> is small, a = 1 — (1 + 
q) 3 4> and 4tt(1 4- q) 3 R 3 /3 is the excluded volume connected 
with each sphere. As <f> increases, the average distance between 
adjacent spheres decreases, and a becomes larger than 1 - (1 + 
q) 3 4> because the depletion layers overlap. 

Because c\ Y /c\ = a(q, <(> n )/a(q f <^ x ), we substitute the ex- 
pression for a. provided by Eq. 5 into the ratio of the free volume 
fractions, and thereby we can determine an apparent q value for 
each experimental tie-line. We find that, within the experimental 
error, this quantity does not depend on the tie-line position in the 
coexistence curve. We present the average value of the apparent 
q as a function of the square root of PEG molecular weight as 
the filled circles in Fig. 3. We note that the q values calculated 
by using a = 1 - (1 + q) 3 <f> are only ~10% smaller than the 
values calculated by using the complete expression for a and so 
conclude that the overlapping of depletion layers contributes 



only marginally to the apparent q values reported in Fig. 3. We 
have also considered the effect of the PEG molecular weight 
polydispersity. Because the polydispersity of PEG is given by a 
narrow Poisson distribution function to be applied to PEG (38, 
39), it can be shown that the polymer distribution in each of the 
two coexisting phases is essentially the same as in the original 
homogeneous system, and hence polydispersity does not signif- 
icantly modify the slope of the tie-lines. 

In Fig. 3, we also present the theoretical q values (open 
circles), calculated by using q = l.\R g /R. R g for PEG is taken 
from ref. 15, whereas R is estimated from the molecular weight 
and specific volume of 7D. As seen in Fig. 3, both theoretical q 
and the apparent q vary as the square root of the PEG molecular 
weight as expected for a Gaussian coil (40). However, the 
apparent values of q are numerically ***50% smaller than the 
corresponding theoretical values. This significant discrepancy 
represents the deviation of the actual 7D-PEG interactions 
relative to the sphere-ideal coil case. 

It has been found that q R g /R in the case of bovine serum 
albumin, chymotrypsinogen, and RNase, whereas q « 0.5R g /R 
in the case of j3-lactoglobulin and lysozyme (12, 15). The 
difference between the values obtained for the various proteins 
suggests that the departure from the excluded volume model is 
related to direct interactions of PEG with protein surface 
residues (15). Our results clearly suggest the presence of such 
specific interactions in the yD-PEG-water systems. It appears 
that weak attractive interactions between PEG and yD induce a 
flattening of the polymer coils near the protein surface and a 
corresponding reduction of the depletion layer thickness. 

The Spinodal Boundaries and the Dependence of the LLPS Tempera- 
ture on PEG Concentration. The spinodal surface (40), which 
describes the spinodal temperature, T sp , as a function of protein 
and PEG concentrations, defines the boundary between the 
stable domain ((d/ti/dci)^ > 0) and the unstable domain 
((d(Li/dCi)fa tT < 0) of a homogeneous protein-PEG-water 
system. The spinodal condition, (3 Ai/dcO^np = 0> can be used 
to determine the spinodal temperature, T sp , as a function of c\ 
and C2. Indeed, on differentiating Eq. 3a with respect to Ci at 
constant fa and T, we obtain: 




In the case of hard sphere-ideal polymer coil excluded volume 
interactions, a = exp(-f), and the difference (d£/dci)r ~ 
(d 2 €/dc\)T becomes equal to the quantity, (d 2 a/dCi) T /cc. If we 
use Eq. 5, we find that, when both and q are small, 
(d 2 a/d</> 2 ) r /a is approximately equal to 12q 3 . Although Eq. 5 
is inaccurate insofar as high-order derivatives such as 
(d 2 a/d<f> 2 )r are concerned, it usefully describes the general 
features of the excluded volume model. Indeed, due to the 
overlapping of adjacent depletion layers, (d 2 oc/d<j> 2 ) T must be 
positive and must increase as the depletion layer thickness 
increases. These general features of the excluded volume model 
can be used to understand the effect of PEG on initially stable 
protein aqueous solutions. In fact, because (df/dci) 2 - — 
(d 2 g/dcf)T is positive, Eq. 6 predicts that as the polymer 
concentration increases, (d(*>i/dCi)^ T decreases and can be 
made to reach zero, i.e., a spinodal boundary. Thus, sufficient 
PEG can in principle induce phase separation at any given 
temperature. This is true even if the protein-protein interactions 
are repulsive. 

We now use the spinodal condition to analyze the dependence 
of the LLPS temperature on PEG concentration. We choose the 
spinodal condition because the spinodal surface is tangent to the 
coexistence surface on the critical line, and because it is math- 
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ematically more convenient than the coexistence condition. If we 
apply the spinodal condition to Eq. 6 and differentiate the 
second member with respect to both C2 and T sp at constant c\, 
we obtain in the limit of C2 approaching zero: 

dc 2 ) Ci ~ (dtfo/dcjTv) ' [?1 

The quantity (d 2 {i , i/dc 1 dT sp ) is equal to ~(d 2 e' /dcl) Tsp /R7%, 
where e' is the internal energy per unit volume of the protein- 
water binary system. On addition of another protein molecule to 
the system, the mean distance between the proteins must 
decrease, promoting more protein-protein interactions. From 
this second-order effect, it follows that the quantity (3 2 « '/<te})r 
is negative for attractive protein-protein interactions, zero for 
hard-core interactions (i.e., hard spheres), and positive for 
repulsive interactions. 

Eq. 7 correctly describes the sign of the dependence of the 
LLPS temperature as a function of PEG concentration for the 
yD-PEG-water systems. The LLPS temperature increases as 
the PEG concentration increases because of the attractive nature 
of the protein-protein interactions [(d 2 £i/ dcidT sp ) > 0]. Also, 
(dT sp /dC2)ci increases with PEG molecular weight because 
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(d£/dci)r - (d 2 £/dc 2 ) T increases as the depletion layer thick- 
ness increases. 

Conclusion 

Our experimental findings and theoretical analysis demonstrate 
that liquid-liquid phase transition in ternary protein-PEG- 
water systems is a powerful tool for understanding the effect of 
PEG as a precipitating agent. We have shown that such studies 
can be used to locate the phase boundary for binary protein- 
water systems, which is very useful if this phase separation is 
hidden. We can also predict the change of protein solubility as 
a function of the concentration of crystallizing agents. Moreover, 
the slope of the tie-lines and the dependence of LLPS temper- 
ature on polymer concentration provide a sensitive check of the 
validity of excluded volume models for protein-polymer inter- 
actions. The increase of the LLPS temperature with PEG 
concentration is due to attractive protein-protein interactions. 
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